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M. Mori 1 , N. Afzal Shooshtary 1 , T. Tohyama 2 , S. Maekawa 1 ' 3 

1 Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
2 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
3 JST, CREST, 5 Sanbancho, Chiyoda-ku, Tokyo, 102-0075, Japan 

The 63 Cu nuclear quadrupole resonance (NQR) frequency, vq, in the multi-layered 
cuprates is calculated in the cluster model by the exact diagonalization method. The charge 
imbalance between the outer CuC>2 plane (OP) with apical oxygen Oa and the inner plane 
(IP) without Oa in three-layered Tl 2 Ba2Ca 2 Cu30io is estimated by comparing our results 
with the experimental vq. In Tl-based cuprates with more than three layers, we predict a 
large enhancement of the splitting of vq between OP and IP by taking into account the 
reduction of bond length between Cu and Oa and a resulting enhancement of the charge 
transfer energy. Our results show that the NQR frequency is a useful quantity to estimate 
the charge imbalance in the multi-layered cuprates. 

KEYWORDS: nuclear quadrupole resonance, multi-layered cuprates, exact diagonalization, api- 
cal site, Madelung potential 

Much study has aimed the higher superconducting transition temperature, T c . Regarding 
the high-T c cuprates, maximum T c depends on materials, e.g., T c ~40 K for La-family and 
T c ~90 K for Y-family, even if carrier densities are optimized. In particular, among hole doped 
cuprates, such a material dependence closely correlates with an energy difference between 
in-plane and apical oxygens, Aea- 1 ^ Actually, a spatial variation of the pairing gap due to a 
modulation of the apical oxygen, Oa, is observed by the scanning tunneling microscopy. 2 ' 3 ) 
In addition, it is known that disorder around the apical site is much more harmful on T c than 
other randomness. 4 ) 

On the other hand, T c of multi-layered cuprates correlates with the number of Cu02 planes 
in a unit cell, n, and shows a maximum at n=3. 5 ' 6 ) One of the reasons why T c is suppressed 
for n >3 is the charge imbalance, i.e., a hole concentration in the pyramidally-coordinated- 
outer-planes (OP's) is different from that in the square-coordinated-inner-planes (IP's). 5 ' 7-10 ) 
The nuclear magnetic resonance (NMR) study has reported that the hole concentration in 
OP, <5op, is larger than that in IP, <5rp. 10 ) The magnitude of charge imbalance, A5= 5op — <5ip, 
increases with n, and becomes 0.10 at most for n=3. 10 ) Such a charge imbalance induces some 
interesting phases, e.g., two kinds of superconducting gap 9 ) and coexistence of superconducting 
and anti- ferromagnetic states. n ~ 16 ) 

To estimate the local charge density, the nuclear quadrupole resonance (NQR) frequency, 
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z/q, is a useful quantity, 17 ' 18 ) since fQ is determined by a local charge distribution in the ground 
state. In the NQR measurement for a three-layered Tl2Ba2Ca2Cu30io (T12223), 18 ** the charge 
imbalance has been observed as two separated peaks with about 6.6 MHz difference. Here, we 
note that the Knight shift, K c , is also proportional to a local charge density on a measured 
site. However, one needs to examine not only the ground state but also the excited states and 
the neighboring spin-spin interactions, since K c is proportional to the susceptibility and the 
transferred hyperfine couplings. 19 ) In addition, one needs to measure not only 63 Cu site but 
also O site to estimate the hole densities in Cu02 plane. However, the NMR measurement 
on 17 site is known to be rather difficult. Hence, an empirical relation between K c on 63 Cu 
at room temperature and the hole density estimated by the NQR measurement has been used 
in general. 10 ) 

Previous theoretical studies on the NQR showed material dependence of vq by neglecting 
a contribution from Oa- 20-22 -* By following Ohta et al. 22 ^ we tried to estimate the charge 
imbalance observed in the NQR measurement by assuming Z\<5=0.10 and no difference of 
charge transfer energy, A, between IP and OP. However, our estimated value of the difference 
was only 2.8 MHz at most. The experimentally observed 6.6 MHz difference is hard to be 
explained without taking account of Oa- 23 ^ To obtain more accurate value of charge imbalance, 
which correlates with T c , one cannot ignore Oa- 

In this Letter, z/q on a Cu site is numerically studied in a cluster model including the 2p z 
orbital of Oa by the exact diagonalization method. We carefully study effects of Aea and A 
on z/q, and show the doping dependence of z^q in IP and OP. To make a link between our 
theoretical result and real multi-layered material, we calculate the Madelung potential of ionic 
model for several different series of multi-layered cuprates. We predict a large enhancement 
of the splitting of isq between OP and IP in Tl-based cuprates with more than three layers. 

The hole densities on the Cu site, which determine z/q, are numerically calculated by 
the exact diagonalization method in the cluster model as shown in Fig. 1(a) for OP. On 
the other hand, for IP, we adopt another cluster, in which 2p z and 4p z orbitals are omitted 
from the cluster in Fig. 1(a). 22 ) We assume that the on-site Coulomb interactions [7^=8.5 
eV, U p =U Pz =4.1 eV and 17^=17^^2.0 eV are material independent, 22 ' 24 ) and the hopping 
integrals 2^=1.06 eV, T pp =0.46 eV, T PPz =0.22 eV, 2^=1.07 eV and 2^^=0.38 eV are 
adjusted to the case of T12223 (ref. 24) by the standard bond-length dependence of T p & cx d~ A , 
Tpp oc d~ 3 with the parameters given in refs. 1 and 21. The notations, p, p z , 4p, 4p z and d, 
denote 2p X;V , 2p z , 4p x>y , 4p z and d x 2_ y 2, respectively. Figure 1(b) is the schematic picture 
of level separations; Aea between p z and p, A between p and d, and A' between d and 4p. 
We assume that 4p x>y and 4p z are degenerate. For both IP and OP, A' is fixed to 4.8 eV. 22 ) 
Regarding the doping rate, S p (p=OP, IP), we calculate two cases of hole number in the 
cluster: 7 and 8 holes (5 and 6 holes) in OP (IP), where 7 (5) holes are associated with the 
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Fig. 1. Schematic pictures of (a) CuOs cluster with Cu Ap a orbitals for OP, and (b) its energy level 
scheme for holes. For IP, CuCU cluster with Cu Ap x { y ) orbitals is adopted, and then 2p z and Ap z 
orbitals are omitted. 22 ) 



half-filling, <5op=0 (<5rp=0). The case of 8 (6) holes in OP (IP) corresponds to the hole-doped 
case with 5op=1 (<5ip=1)- Since the cluster size is small, we assume that an intermediate case 
of 5 P could be linearly interpolated between z/q's of <5 P =0 and 1: 

v Q (5 p ) = (1 - S p )vq(6 p = 0) + v Q (5 p = l)5p, (1) 

where the notations, z/q(#ip) and z/q(5op)> are used to denote the NQR frequency in IP with 
5\p and OP with <5op, respectively. 

On a 63 Cu nucleus with the quadrupole moment Q = —0.211 barn and the spin angular 
moment / = 3/2, vq in OP is calculated by 

eQ ( 

= 2^'\ n 3cl 1( a_ v a Af 3d + (4 - n 4px - n Apy )M 4p 

- 2(2 -n 4pz )M 4p y (2) 

The hole density on a Cu site is denoted by n« (i = 2>d x 2_ y 2,Ap x ,Ap y , 4p 2 ). Since the third term 
in eq. (2) represents the Oa contribution to OP, vq in IP is given by the sum of the first and 
the second terms. We use the values M 3d =-47 x 10 21 V/m 2 and M 4p =125xl0 21 V/m 2 . 22 ) 
The M4 p -terms in eq. (2) arise from Ap a (a=x, y, and z) orbitals on a Cu site, which are 
partial- wave components of the 2p a wave function in the neighboring O sites. 21 ' 22 ) The z-axis 
is perpendicular to the plane, and the isotropy of x and y is assumed. 

In Fig. 2, z/q is plotted as a function of A for a fixed value of Aea=2.5 eV. We note that 
moderate variation of Ae\ does not change the results so much. The symbols shown in Fig. 2 
are fitted by a quadratic form, vq = Aq + A±A + A2A 2 . We find that the values of A\ and A2 
between OP and IP are similar for a given <5 p . 25 ) vq increases monotonically with A, because 
of the increase of n^d 2 2 as reported in Ref. 21. The large enhancement of vq (~ 25 MHz) 
from 5 P =0 to 5 P =1 is also attributed to the increase of n^d a 2 - -^ or a gi yen 8 P , z/q at OP 
is larger by ~ 2.5 MHz than z/q at IP. This is due to the contribution of the third term in 
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Fig. 2. The calculated NQR frequency vq as a function of A. Square and triangle symbols are the 
results for IP and OP, and the open and closed symbols mean the undoped ((5 p =0) and doped 
(<5 p =l) cases, respectively. These symbols are fitted by a quadratic form uq = Aq + A\A + A 2 A 2 . 



eq. (2). 

It is worth emphasizing that Fig. 2 is available for evaluating one of the three quantities 
(Vq, A, and 5 P (see eq. (1)) when two of them are provided. As an example, let us evaluate the 
charge imbalance A5 for T12223 to which the 6.6 MHz difference of vq's between OP and IP 
has been reported. First of all, we need to know A at OP and IP. Since it is very difficult to 
evaluate A experimentally, we have to rely on theoretical results. As will be discussed below, 
both A's at OP and IP are estimated to be ~2.5 eV when n < 3. Putting the estimated A 
into the fitted equation of vq, we find 

vq(Sop) - fq(<5rp) ~ 3.0 + 24.0A5. (3) 

This is a useful relation between the frequency difference and charge imbalance for the case of 
n < 3. The 6.6 MHz difference of vq gives A<5=0.15, which is slightly larger than 0.10 estimated 
by the empirical relation used in the NMR study. 10 - 1 One possibility of this discrepancy is that 
the linear interpolation eq. (1) is too simple. Another one is that the empirical relation may 
need to be modified by considering the coordination number of oxygens around Cu site. 

We can see in Fig. 2 that vq sensitively changes with A, i.e., 0.1 eV increment of A provides 
5~6 MHz enhancement in isq. Therefore, to study the charge imbalance in the multi-layered 
cuprates, we need to estimate A in each plane. Below, the following three series of multi-layered 
cuprates are examined (1) Tl-12 series, TlBa 2 Ca n _iCu n 2n + 3 , 26 - 29 ) (T C =48K, 110K, 123K, 
112K, 107K from n=l to n=5), (2) Tl-22-series, Tl 2 Ba 2 Ca n _iCu n 2n+4 , 29 - 33 ) (T C =95K, 
118K, 125K, 112K, 105K from n=l to n=5), (3) Hg-12-series, HgBa 2 Ca„_iCu„0 2ri+2 , 34 - 37 ) 
(T C =98K, 128K, 135K, 127K, 111K from n=l to n=5). The distance between Cu and A , d A , 
is plotted as a function of n in Fig. 3. Obviously, d A decreases for n > 3 except for Hg- family. 
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Fig. 3. The bond length between Cu and Oa, <^a, versus the number of CuC>2 layer, n, for three series 
of multi-layered cuprates. The atomic coordinates are found in the following references: Refs. 26 ~ 29 ) 
for Tl-12 family, Refs. 29 " 33 ' for Tl-22 family, Refs. 34 - 37 ) for Hg-12 family. 




Fig. 4. Potential difference between apical oxygen and in-plane oxygen, Aea= (V(Oa) — V(0))/eo, 
versus n for three different series of multi-layered cuprates. 



This may be a common feature of Tl-based cuprates. 38 ) The suppression of d\ may be relevant 
to the suppression of T c , since the apical oxygen position directly correlates with Aea and the 
paring gap. 1,3 ) 

To estimate As a and A, we calculate the Madelung potential V(A) (A=Cu, Oa, and O 
in the plane) in the ionic model. In Fig. 4, Aea= (V(Oa)— V(0))/eo is plotted with n for 
each family. Note that the screening effect by the long-ranged Coulomb interaction can be 
effectively expressed by dividing the Madelung potentials by a constant eo, whose value is 
assumed to be eo = 3.5. 1 ' 39 ' On the whole, Aea does not strongly depend on n, while each 
family shows different trend. In Tl-12 families, As a slightly decreases with n (> 2), while 
Aea in Hg-family increases up to n = 4. Since the moderate variation of Aea shown in Fig. 
4 does not change uq so much, we can assume Aea=2-5 eV in calculating z/q. 

On the other hand, one cannot ignore the variation of A shown in Fig. 5, where A in OP, 
A(OP), is plotted as a function of n. We used the relation, A = (V(O) - V(Cu))/e - 10.9. 1} 
It is noticed that A(OP) of Tl-families increases particularly for n=4 and 5. The deviation of 
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Fig. 5. The charge transfer energy in the outer plane, A(OP), versus n for three series of multi- 
layered cuprates. The charge transfer energy in the inner plane, A(IP), is almost independent of 
n as, A(IP)~2.5 eV. 

Tl-family or the flatness of Hg-family must be due to their characteristic local structure as 
shown in Fig. 3. Note that A in IP, A(IP), is always about 2.5 eV independent of family. Let 
us consider T11245, i.e., n=5 in the Tl-12 family. Taking A(OP)=3.4 eV and A(IP)=2.5 eV, 
we can predict from Fig. 2 a large NQR splitting, i/q(OP)— z/q(IP)~25 MHz, for the case 
of no carrier (<5op=^ip=0). Since real T11245 has carriers and <5op>^iP; the NQR splitting 
for T11245 must be more enhanced than 25 MHz. Such a large NQR splitting has not been 
observed so far, but careful NQR measurements may discover a large enhancement of the NQR 
splitting in the Tl-families with n=4 and 5. When we try to estimate the charge imbalance 
from observed fQ, we should be careful about the difference of A(OP) and A(IP). For example, 
for T12234 or T12245 one should take A(OP)w 3.0 eV and A(IP)« 2.5 eV according to Fig. 5 
in obtaining a relation between the splitting and the charge imbalance. The relation reads 
^q(^op) - ^q(^ip) = 17.64 + 24.0(5 O P - ftp) - 1.5«5 p ~ 17.64 + 24.0A<5. If one observed a 
NQR splitting, for instance, with 20 MHz, one might obtain A<5=0.10. 

To conclude, we studied the NQR frequency, vq, in the cluster model including the con- 
tribution from p z orbital of apical oxygen by the exact diagonalization method. The charge 
imbalance in the three-layered T12223 is estimated by comparing our result and the experi- 
mental uq. Our estimate on the charge imbalance A5=5op— 5ip is 0.15, which is larger than 
the NMR result, A5=0.10. The Madelung potential in several types of multi-layered cuprates 
is calculated in each plane to show the n-dependences of energy difference between apical and 
in-plane oxygens, Aea, and charge transfer energy, A. It is found that a large splitting of 
fQ between OP and IP can be expected in the Tl-12 and Tl-22 families. 

Recently, the angle-resolved photoemission spectroscopy (ARPES) study in 
Bi2Ba2Ca2Cu30io has reported that <5(OP)=0.26 and 5(IP)=0.06, which is a quite 
large charge imbalance, Z\<5=0.20. 40 ) In the ARPES study, the hole densities are estimated by 
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the ratio of the Fermi surface volume and the first Brillouin zone. The ARPES result is larger 
than both ours and the NMR result. These discrepancies among theory and experiments 
must be solved in near future. 
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